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Constructing Parton Convolution in Effective Field Theory 
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Parton convolution models have been used extensively in describing the sea quarks in the nucleon 
and explaining quark distributions in nuclei (the EMC effect). From effective field theory point of 
view, we construct the parton convolution formalism which has been the underlying conception of 
all convolution models. We explain the significance of scheme and scale dependence of the auxiliary 
quantities such as the pion distributions in a nucleon. As an application, we calculate the complete 
leading nonanalytic chiral contribution to the isovector component of the nucleon sea. 



Calculating parton distributions from the first princi- 
ples have been proven difficult. The only approach at 
present is lattice field theory in which the moments of 
parton distributions are simulated on a Euclidean lat- 
tice pl-H. A useful phenomenological approach used to 
understand certain aspects of the parton distributions, 
such as their modifications in a nucleus or the origin of 
sea quarks, is to construct a parton convolution model 
|^|~H. The essence of the parton convolution picture is 
easy to describe: First, the system under consideration is 
treated as a composite of hadrons (nucleons and their ex- 
citations and mesons). Then hard scattering is pictured 
to happen on one of its hadron constituents in which 
the parton distributions are presumed known. Although 
widely used, the convolution models seem hard to justify 
at a more fundamental level. For instance, it has been 
difficult to construct systematic corrections to the model 
predictions. 

In recent publications p,pi3[, the effective field theory 
techniques — more specifically chiral perturbation theory 
| PT| , P^ — were first used to understand the chiral correc- 
tions to the parton distributions in the nucleon. We find 
that the approach can easily be generalized to under- 
stand the sea quark distributions in the nucleon and the 
modifications of parton distributions in nuclei. The re- 
sult is a rigorous and general way to construct parton 
distributions of composite systems out of their hadron 
constituents, generalizing the conventional convolution 
model approach. Unlike the convolution models, how- 
ever, effective field theories allow a systematic way to 
account for higher-order contributions through power 
counting. 

Our observation is simple: The moments of the parton 
distributions are defined from the matrix elements of the 
twist-two operators. For instance, for the unpolarized 
parton distributions these operators are 






97 



(''HDi'^ ■■■iD'^^h 



(1) 



In effective field theories, these operators are matched to 
hadronic operators with the same quantum numbers. 









where j labels different types of hadronic operators and 
Cg" are c- number coefficients outside of the effective field 
theory framework. As a specific example, the isovector 
operator 
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can be matched onto the hadronic operators 



O^ 



aj-^ JTr SVgiD^^i • • • iD^'^^T. 
+Tr T.TsiD'-t'^ •••iD^")Et | 
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where N and N are the nucleon fields in the nonlin- 
ear representation, v^ is the nucleon four-velocity, u = 
exp{iTr°'T"' /2fTr) with pion fields 7r° and decay constant 
/tt — 93 MeV, S = u^, and (• • •) denotes the symmetriza- 
tion of the indices in between. We have neglected terms 
with more derivatives as well as with multiple nucleon 
fields relevant for the nuclear modifications of parton 
distributions. 

Of course, the expressions in Eqs. (g) and (^ are use- 
less if there are no small expansion parameters and every 
term contributes at equal importance. If, however, a sen- 
sible expansion scheme exists, as when an effective field 
theory becomes applicable, one can learn some aspects 
of the parton distributions through the expansion. In 
particular, the expansion allows the quark distributions 
in a hadron system H be expressed as convolutions; 



Qnix) = Y^ 
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(2) where qj{y) is defined through its moments 
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and can be interpreted as the quark distribution in a 
hadron state j. fjHiv) is defined through moments 



rfyy"-V,ff(y) = 



1 



2{p+y 



-(P|0+-+|P)h , (7) 



which can be interpreted as the hadron j distribution 
in H (P+ = {P° + P^)/y/2, and similarly for other in- 
dices). Using Eqs. (||)-(0), it is easy to check that the 
(n— l)th moment of qnix) reproduces Eq. (0). A similar 
convolution formula can be derived for polarized distri- 
butions. Equation (||) is schematic because j refers to 
general hadronic states having the quantum numbers of 
the twist-two operators. For example, j can be an iV — A 
interference state, or a 3-pion state. For intermediate 
states with multiple hadrons, y must be extended to sev- 
eral light-cone variables. For off-diagonal hadron states, 
the quark distributions are the generalized parton distri- 
butions that have been studied recently in the literature 
|l3| . According to Eq. (||), various convolution models 
consist of particular truncations of the expansion. 

To illustrate the above formalism, we study in the re- 
mainder of the paper the moments of the proton's isovec- 
tor anti-quark distribution, up{x) — dp{x), in chiral per- 
turbation theory. A chiral contribution to this quantity 
was first computed in Ref. ||lj|, and our result partially 
confirms the answer there. However, we have also found 
an additional contribution which depends on the quark 
distribution in the chiral limit. 

We start by expressing the isovector, spin- averaged 
twist-two quark operators in terms of the operators with 
pure fields, 
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where after the second equal sign, we have made the 
chiral expansion and neglected the higher-order terms 
in r7T,7r/(47r/7r). The matrix elements of the isovec- 
tor operators, O'^X-'^'- = —ie^^-^TT^id'^^---id'^"TT^, in 
the proton state define the isovector pion distribution 

fi7T+-7T-)/p{y)'- 
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The matrix elements A„ can be calculated in chiral 
perturbation theory. In dimensional regularization and 
heavy- nucleon formalism [|ll| , we find 
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where gA is the neutron decay constant in the chiral 
limit. The explicit nucleon mass dependence comes from 
the kinematic factors in the definition of the matrix el- 
ements and from setting the renormalization scale (in 
the chiral logarithm) to mjv. The ellipses denote terms 
analytic in quark masses or subleading in chiral power 
counting. The n = 1 moment is special because it re- 
ceives a pion tadpole contribution. The Bose symmetry 
leads to vanishing n = even moments. 

From the moments, one can construct the pion distri- 
bution, /(T+-ir-)/p(2/) in the proton. 
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where 
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Apart from the leading nonanalytic contribution from 
small y, the above expression also contains analytic 
terms in quark masses. The Bose symmetry dictates 
f{Tr+-Tr-)/p{y) as an even function of y. Equation (O) 
can also be derived from the definition of the pion dis- 
tribution as the matrix elements of the light-cone string 
operator in the same scheme p5[ 
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dX 
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where n ~ (1,0,0,-1) is a light-cone vector. The 5- 
function contribution at y = comes from the pion- 
nucleon seagull vertex. 

Before going further, it is useful to discuss the scheme 
and scale dependence of the pion distribution in the nu- 
cleon. For simplicity, we match a generic quark operator 
Oq schematically onto the sum of a pion operator O^ 



and a nucleon operator O 



N, 



Ogifi) = a{fi, A)0,(A) + 6(/i, A)Ojv(A) 



(14) 



The quark operator must be renormalized in QCD and 
has a perturbative QCD renormalization scale /i much 
larger than the QCD scale, Aqcd- The fi dependence 
in the right-hand side of the equation is confined to the 
coefficient functions a(/i,A) and b{fi,A). The hadronic 
operators, on the other hand, have an implicit scheme 



and an explicit renormalization scale A (< 1 GeV) de- 
pendence in chiral perturbation theory. The scheme and 
scale dependence in these operators ought be cancelled 
entirely by the scheme and scale dependence of the co- 
efficient functions. When Eq. (nj) is sandwiched in the 
nucleon state, (N\OTr\N) is interpreted as the pionic ef- 
fect in the nucleon. Because of the scheme and scale 
dependence of the matrix element, we cannot define in 
an absolute sense, for instance, the number of pions in 
the nucleon. 

In a typical calculation of the Sullivan process, the 
baryons (nucleons) are treated fully relativistically with 
finite masses (see, for example, Ref. |14 ). The leading- 



order contribution is proportional to the coupling 
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where we have made use of the Goldberger-Treiman re- 
lation. Since the nucleon mass toat is on the order of 
47r/^, the expansion parameter is by no means small. 
In particular, the multi-loop contributions are propor- 
tional to (5,rA'Af/47r)^" which is not suppressed relative 
to the leading order. Therefore, in this scheme of chiral 
expansion, {N\OTr\N) cannot be calculated with a con- 
trolled approximation. As such, the leading order result 
can only be taken as a model prediction even when the 
loop momentum of the pion is kept small. The heavy- 
baryon chiral perturbation pl| is designed to overcome 
this difficulty. In this scheme, the leading-order result is 
quadratically divergent. In a cut-off regularization, it is 
proportional to g^A"^ /{ATrf-^}'^, which is a small param- 
eter when A <C 47r/^. As a consequence, the pionic con- 
tribution is very sensitive to the cut-off. In dimensional 
regularization, on the other hand, the power divergences 
are absent by definition. From dimensional analysis, the 
matrix element (iV|0^|iV), and here the pion distribu- 
tions in the nucleon, vanishes in the chiral limit. The 
renormalization scale A enters in the logarithms of the 
leading nonanalytical term and controls the relative con- 
tribution to Oq from O^ and On- 

The coefficient a„ in Eqs. (Q) and (||) is related to the 
u — d distribution in the 7r+ meson: 
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where the superscript indicates the chiral limit. The 
isospin charge counting gives ai — 2. The difference 
between u*^ , and the real- world u«.+ is a chiral correction 
which has been computed in Ref. |Ifl]: 
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where the delta function contribution comes from isospin 
conservation. For the following purpose, this correction 
is of higher order. 

Through Eq. (H), we construct a contribution to the 
proton's u{x) — d(x) distribution from the pion cloud: 

u^/p{x)-d^/p{x)= / —\ul+{y)-dl+{y) 

J\x\ V ^ 



-(4+(2/)-d^+(y))}/(^+_^-)/, 
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where we have used the symmetry property of 
f(Tr+-TT-)/N{y) in y and the definition of antiquark dis- 
tribution qj{y) — ~qj{—y). Since the contribution is an 
even function of x, we can easily calculate the resulting 
up-down asymmetry in the sea. 
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Above result agrees with that in Ref. [Q in the limit 
gA ~ i- A number of comments can be made about this 
result. First, L dx(d{x) — u{x)) is not an usual moment 
of the quark distribution and cannot be calculated as the 
matrix element of a local operator in QCD or effective 
field theories. We overcome this difficulty by converting 
the matrix elements of the complete tower of twist-two 
operators into a distribution and then integrating over 
X. The simple result follows from that the pion contri- 
bution to the quark distribution in the nucleon is even 
in X. Therefore, in a certain sense, our derivation pro- 
vides a formal justification for the approach used in Ref. 
[H. Second, the difference between our result and that 
of Ref. iQ comes from the use of the linear sigma model 
in the later paper. As was discussed in Ref. 0], in a gen- 
eral formulation of chiral expansion in which the nucleon 
fields furnish a linear representation of the chiral group, 
the difference disappears. Finally, the above nonanalytic 
chiral contribution to /p dx{d(x) — u{x)) has an oppo- 
site sign from the experimental data ||l^ , indicating the 
importance of other contributions to the nucleon sea. 

According to the expansion in Eq. (Q), the twist-two 
pion operators are not the only source of chiral log- 
arithms in the antiquark distributions in the nucleon. 
They can also be generated from the nucleon operators. 
Using vPp{x) and dp (a:) to represent the quark distribu- 
tions in the proton in the chiral limit, the result in Refs. 
[p|,pi)[ allows us to find another chiral contribution to u 
and d distributions in the proton, 

un/p{x) - dN/p{x) 
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Integrating from —1 to 0, we have 
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The second term in the bracket represents a chiral con- 
tribution which cannot be calculated without the knowl- 
edge of Jn{d%{x) — u°p{x))dx. Together with the result 
in Eq. (p8|), we have the complete leading chiral con- 
tribution to the isovector antiquark distributions in the 
nuclcon. 

The leading chiral contribution to the polarized up 
and down quark asymmetry in the proton sea can be 
calculated in a similar way. The quark-helicity up-down 
asymmetry in the sea is 
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where the superscript indicates the chiral limit. Simi- 
larly, the quark-transversity up-down asymmetry in the 
sea is 



{Sdp{x) — Sup{x))dx — / {6d!j,{x) — 5u^p{x))dx 
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In both cases, the contribution from the pion cloud is of 
higher order in chiral power counting. 
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